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We study the growth of perturbations in an expanding Newtonian universe with Bose-Einstein condensate dark matter. 
We first ignore special relativistic effects and derive a differential equation governing the evolution of the density 
contrast in the linear regime taking into account quantum pressure and self-interaction. This equation can be solved 
analytically in several cases. We argue that an attractive self-interaction can enhance the Jeans instability and fasten the 
formation of structures. Then, we take into account pressure effects (coming from special relativity) in the evolution 
of the cosmic fluid and add the contribution of radiation, baryons and dark energy (cosmological constant). For a 
BEC dark matter with repulsive self-interaction (positive pressure) the scale factor increases more rapidly than in the 
standard ACDM model where dark matter is pressureless while for a BEC dark matter with attractive self-interaction 
(negative pressure) it increases less rapidly. We study the linear development of the perturbations in these two cases 
and show that the perturbations grow faster in a BEC dark matter than in a pressureless dark matter. This confirms a 
recent result of Harko (2011). Finally, we consider a "dark fluid" with a generalized equation of state p = (ap + kp 2 )c 2 
having a component p — kp 2 c 2 similar to a BEC dark matter and a component p = ape 2 mimicking the effect of the 
cosmological constant (dark energy). We find optimal parameters that give a good agreement with the standard ACDM 
model assuming a finite cosmological constant. 
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1. Introduction 

Several astrophysical observations of distant type la super- 
novae have revealed that the content of the universe is made 
of about 70% of dark energy, 25% of dark matter and 5% 
of baryonic (visible) matter (Riess et al. 1998, Perlmuttcr 
et al. 1999, de Bernardis et al. 2000, Hanany et al. 2000). 
Thus, the overwhelming preponderance of matter and en- 
ergy in the universe is believed to be dark i.e. unobservable 
by telescopes. The dark energy is responsible for the accel- 
erated expansion of the universe. Its origin is mysterious 
and presumably related to the cosmological constant. Dark 
energy is usually interpreted as a vacuum energy and it be- 
haves like a fluid with negative pressure. Dark matter also 
is mysterious. The suggestion that dark matter may con- 
stitute a large part of the universe was raised by Zwicky in 
1937. He realized that some mass was "missing" in order to 
account for observations. This missing mass problem was 
confirmed later by more accurate measurements (Borriello 
& Salucci 2001). The rotation curves of neutral hydrogen 
clouds in spiral galaxies measured from the Doppler effect 
are found to be roughly flat with a typical rotational ve- 
locity Voo ~ 200km/s up to the maximum observed ra- 
dius of about 50 kpc. This mass profile is much more ex- 
tended than the distribution of starlight which typically 
converges within ~ f kpc. This implies that galaxies are 
surrounded by an extended halo of dark matter whose mass 



M(r) 



rv: 



3 /G increases linearly with radius. Although 



some authors like Milgrom (1983) propose a modification 



of Newton's law (MOND theory) to explain the rotation 
curves of spiral galaxies without invoking dark matter, the 
dark matter hypothesis is favored by most astrophysicists. 

The nature of dark matter (DM) is one of the most im- 
portant puzzles in modern physics and cosmology. A wide 
"zoology" of exotic particles that could form dark matter 
has been proposed. In particular, many grand unified the- 
ories in particle physics predict the existence of various ex- 
otic bosons (e.g. axions, scalar neutrinos, ncutralinos) that 
should be present in considerable abundance in the uni- 
verse and comprise (part of) the cosmological missing mass 
(Primack et al. 1998, Overduin & Wesson 2004). Even if 
the bosonic particles have never been detected in accelera- 
tor experiments, they are considered as leading candidates 
of dark matter and might play a significant role in the evo- 
lution and in the structure of the universe. 

If dark matter is made of bosons, they should have 
formed compact gravitating Bose-Einstein condensates 
(BEC) such as boson stars. Boson stars were introduced 
by Kaup (f 968) and Ruffini & Bonazzola (1969) in the six- 
ties. Early works on boson stars (Thirring 1983, Breit et al. 
1984, Takasugi & Yoshimura 1984, van der Bij & Gleiser 
1987) were motivated by the axion field that was proposed 
as a possible solution to the strong CP problem in QCD. 
For particles with mass m ~ f GeV/c 2 , the maximum mass 
of a boson star, called the Kaup mass, is much smaller 
than the solar mass {Mxaup ~ 10~ 19 A/ Q !) so that these 
mini boson stars, like axion black holes, are not very as- 
trophysically relevant. Some authors (Baldeschi et al. 1983, 
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Sin 1994, Hu et al. 2000) have proposed that dark mat- 
ter halos could be giant systems of "Bose liquid" but in 
that case the mass of the bosons must be extremely small 
(m ~ 10 _24 eV) to yield masses consistent with the mass 
of galactic halos. Such an ultralight scalar field was called 
"fuzzy cold dark matter" (FCDM) by Hu et al. (2000) who 
discussed its overall cosmological behavior. On the other 
hand, Colpi et al. (1986) have shown that if the bosons 
have a self-interaction, then the mass of the boson stars 
can considerably increase, even for a small self-interaction. 
For m ~ 1 GeV/c 2 it becomes of the order of the solar mass 
so that dark matter could be made of numerous boson stars. 
On the other hand, for m ~ leV/c 2 , the mass of the bo- 
son stars becomes of the order of the mass of the galactic 
halos. Therefore, some authors (Lee & Koh 1996, Peebles 
2000, Goodman 2000, Arbey et al. 2003, Bohmer & Harko 
2007) proposed that dark matter halos themselves could be 
in the form of gigantic self-gravitating Bose-Einstein con- 
densates with short-range interactions described by a single 
wave function ip(r,t). In the Newtonian limit, which is rel- 
evant at the galactic scale, the evolution of this wave func- 
tion is governed by the Gross-Pitaevskii-Poisson (GPP) sys- 
tem (Bohmer & Harko 2007). Using the Madelung (1927) 
transformation, the GP equation turns out to be equivalent 
to hydrodynamic (Euler) equations involving an isotropic 
classical pressure due to short-range interactions (scatter- 
ing) and an anisotropic quantum pressure (or a quantum 
potential) arising from the Hciscnberg uncertainty princi- 
ple. For a standard BEC with quartic self-interaction, the 
equation of state is that of a polytrope with index n = 1. 
At large scales, scattering and quantum effects are negligi- 
ble and one recovers the classical hydrodynamic equations 
of cold dark matter (CDM) models which are remarkably 
successful in explaining the large-scale structure of the uni- 
verse. At small scales, gravitational collapse is prevented 
by the (repulsive) scattering or by the uncertainty princi- 
ple. This may be a way to solve the problems of the CDM 
model such as the cusp problem and the missing satellite 
problem (Hu ct al. 2000). 

In our previous papers (Chavanis 2011, Chavanis & 
Dclfini 2011), we performed an exhaustive study of the 
equilibrium configurations of a Newtonian self-gravitating 
BEC with short-range interactions. For a given value of 
the scattering length, we obtained the mass-radius relation 
AI(R) connecting the non-interacting case (corresponding 
to small masses) studied by Ruffini & Bonazzola (1969) 
to the Thomas-Fermi (TF) limit (corresponding to large 
masses) investigated by Bohmer & Harko (2007). We also 
considered the case of attractive self-interaction. This cor- 
responds to a negative scattering length (a s < 0) yielding 
a negative pressure. In that case, we found the existence of 
a maximum mass M max = 1.012fi/ \/\a s \Gm above which 
the system becomes unstable. It turns out that this mass 
is ridiculously small (being possibly as small as the Planck 
mass Mp = 2.1810 _8 kg!) meaning that a self-gravitating 
BEC with attractive short-range interactions is extremely 
unstable (except if the mass of the bosons is extraordinarily 
small like in Hu et al. 2000). We proposed that the scatter- 
ing length could be negative in the early universe and that 
it could help for the formation of structure^. Then, the 



1 This is a potentially interesting idea because, in an expand- 
ing universe, the condensation process is often regarded to be too 
slow to account for the formation of structures (Bonnor 1957, 



scattering length could become positive and help to stabi- 
lize the structures against gravitational collapse. Of course, 
the mechanism by which the scattering length changes sign 
remains to be established so that this idea is highly specu- 
lative. We may note that some atoms in terrestrial BEC ex- 
periments are reported to have negative scattering lengths 
(Dalfovo et al. 1999). Therefore, the possibility of negative 
pressure for a BEC can be contemplated. Furthermore, it 
has been experimentally demonstrated that, under certain 
conditions, it is possible to manipulate the sign and value of 
the scattering length (Fedichev et al. 1996). The extension 
of these ideas to cosmology remains an (important) prob- 
lem that we shall not discuss further here. We shall assume 
that the dark matter in the universe is a BEC with repul- 
sive (a s > 0) or attractive (a s < 0) self-interaction and 
theoretically explore the consequences of this hypothesis. 

If dark matter halos are BECs, they have probably 
formed by Jeans instability. The gravitational instability 
of a scalar field equivalent to a BEC was considered by 
Khlopov et al. (1985), Bianchi et al (1990), Hu et al. (2000) 
and Sikivie & Yang (2009). However, these authors started 
their analysis from rclativistic field equations and did not 
take into account the self-interaction of the particles. In our 
previous paper (Chavanis 2011), we studied the Jeans in- 
stability of a self-gravitating BEC with short-range interac- 
tions described by the Gross-Pitaevskii-Poisson system in a 
static universe. We considered both attractive and repulsive 
self-interaction and found that when the scattering length 
is negative the growth rate of the instability increases with 
respect to the case where the scattering length is positive 
or zero. The next step is to study the Jeans instability of 
a BEC in an expanding universe. This is the topic of the 
present paper. 

In the first part of the paper (Sees. HE]), we neglect 
rclativistic effects and use the equations of Newtonian cos- 
mology introduced by Milne (1934) and McCrea & Milne 
(1934). We study the linear development of perturbations 
of a self-gravitating BEC with short-range interactions in 
an expanding Einstcin-dc Sitter (EdS) universe. In the TF 
approximation, the equation of state of a BEC with quar- 
tic self-interaction is that of a polytrope of index 7 = 2 
(n = 1). In the non- interacting case, we find that the BEC 
behaves similarly to a polytrope of index 7 = 5/3 [n = 3/2) 
but with a quantum mechanically modified Jeans length. In 
the two cases, the equation for the density contrast can be 
solved analytically in terms of Bessel functions. 

In the second part of the paper (Sec. [8]), we take spe- 
cial relativistic effects into account and use the equations of 
Newtonian cosmology with pressure introduced by McCrea 
(1951). In that case, the evolution of the cosmic fluid de- 
pends on the equation of state. We consider a BEC dark 
matter with equation of state p = kpc 2 and add the con- 
tribution of radiation, baryons and dark energy (cosmo- 
logical constant). For a BEC dark matter with repulsive 
self-interaction k > (positive pressure) the scale factor 
increases more rapidly than in the standard ACDM model 
where dark matter is prcssureless (k — 0) while for a BEC 
dark matter with attractive self-interaction k < (negative 
pressure) it increases less rapidly. We study the linear de- 
velopment of the perturbations in these two cases and show 



Harrison 1967). An attractive self-interaction could substantially 
enhance the instability. 



2 



P.H. Chavanis: Growth of perturbations in an expanding universe with Bose-Einstein condensate dark matter 



that the perturbations grow faster in a BEC dark matter 
than in a pressureless dark matter. 

In the third part of the paper (Sec. 0, we consider a 
"dark fluid" with a generalized equation of state p = (ap 
kp 2 )c 2 having a component p — kp 2 c 2 similar to a BEC 
dark matter and a component p = ape 2 mimicking the 
effect of the cosmological constant (dark energy) . We find 
optimal parameters that give a good agreement with the 
standard ACDM model. 

While our scries of papers on this subject was in course 
of redaction, Harko (2011) published a paper where he also 
considered the formation of structures in an expanding uni- 
verse made of BEC dark matter. Our contribution is com- 
plementary to Harko's work and confirms his main results. 
In addition, we address the following issues: (i) we study the 
effect of the quantum pressure in the first part of the paper; 
(ii) we use different relativistic hydrodynamic equations to 
model the cosmic fluid in the second part of the paper; (iii) 
we study a generalized equation of state in the third part of 
the paper; (iv) we consider both positive and negative scat- 
tering lengths throughout the paper. The papers of Bohmcr 
& Harko (2007) and Harko (2011) show that a cosmic BEC 
can be an interesting model for the dark matter of the uni- 
verse. Our series of papers, which develop this idea, go in 
the same direction. 

Finally, in a different perspective, Widrow & Kaiser 
(1993) have proposed to describe a classical collisionless 
self-gravitating system by the Schrodingcr-Poisson system. 
In this approach, the constant Ti is not the Planck con- 
stant, but rather an adjustable parameter that controls 
the spatial resolution XdeB through a de Broglie relation 
AdeB = Ti/mv. It is argued that when h — > 0, the Vlasov- 
Poisson system is recovered and that a finite value of fi 
provides a small-scale regularization of the dynamics. In 
that case, the Schrodingcr-Poisson system has nothing to 
do with quantum mechanics since it aims at describing the 
evolution of classical collisionless matter under the influ- 
ence of gravity (in static or expanding universes). Still, on 
a mathematical point of view, these equations are equiva- 
lent to those describing self-gravitating BECs without self- 
interaction. Therefore, the results of the first part of our 
paper can have application in that context, independently 
of quantum mechanics. 



2. The Gross-Pitaevskii-Poisson system 

Following Bohmcr & Harko (2007), we assume that dark 
matter is a Bose-Einstein condensate. A self-gravitating 
BEC with short-range interactions is described by the 
Gross-Pitaevskii-Poisson system 



dip 
dt 



= -2^ 



m($ + h{p))ip, 



A$ = 4:nGNm\ip\ 2 - A, 



(1) 



(2) 



where p = Nm\ip\ 2 is the density, $ the gravitational poten- 
tial and h(p) = J usR(r—r')p(r' , t) dr' an effective potential 
taking into account small-scale interactions. For the sake of 
generality, we have included the cosmological constant A in 
the Poisson equation^. We write the wave function in the 



2 We may note that the cosmological constant A plays the 
same role as a global rotation of the universe if $ is interpreted 



form ip(r,t) = A(r, i)e lS ( I%t )/ ?i where A and S are real, and 
make the Madelung (1927) transformation 



+ p = Nrn\ip\ 2 = NmA 2 



-VS, 



(3) 



where p(r, t) is the density field and u(r, t) the velocity 
field. We note that the flow is irrotational since V x u = 0. 
With this transformation, it can be shown that the Gross- 
Pitaevskii equation ([1]) is equivalent to the barotropic Euler 
equations with an additional term Q called the quantum 
potential (or quantum pressure). Indeed, one obtains the 
set of equations 



dp 
dt 



V ■ (pu) = 0, 



du , , 1 1 
— + u • V)u = — Vp - V$ VQ, 

dt p m 

A$ = 4nGp - A, 



with 



Q = - 



h 2 A^p 

2m y/p 



(4) 



(6) 



(7) 



The pressure p(p) is determined by the effective poten- 
tial h(p), playing the role of an enthalpy (Chavanis 2011), 
through the relation p'(p) = ph'(p). For a contact pair in- 
teraction usn(r — r') = gS(r — r'), the effective potential 
h(p) = gp where g = Ana s h /m 3 is the pseudo-potential 
and a s is the s-scattering length (Dalfovo et al. 1999). For 
the sake of generality, we allow a s to be positive or negative. 
The corresponding equation of state is 



2-KOLeTi 2 



(8) 



It corresponds to a polytropic equation of state of the form 



p = Kp 1 



7 = 1 + -, 



(9) 



with polytropic index n = 1 (i.e. 7 = 2) and polytropic 
constant K = 2-Ka 8 Ti 2 jrr? . 

The effective potential corresponding to the polytropic 
equation of state (|9|) is 



Hp) 



7-1' 

This leads to a GP equation of the form 



ih 



dip 
~dt 



2m 



AiP + m(<i> + K\iP\ 2/n )iP, 



(10) 



(11) 



where k = K(n + \){Nm) x l n . This is the usual form of 
the GP equation considered in the literature (Sulem & 
Sulcm 1999). The standard BEC, with a quartic inter- 
action p (x p 2 cx \ip\ 4 , corresponds to n = 1. We may 
also recall that classical and ultra-relativistic fermion stars 
are equivalent to poly tropes with index n — 3/2 and 
n = 3 (Chandrasekhar 1939). The GPP system (|TTj) -(|2"]) 



as an effective gravitational potential $e// = $ — X r) 2 
in the rotating frame. Indeed, the Poisson equation becomes 
A$ e // = 4tyGp — 2Q 2 which allows the identification A = 2$1 2 . 
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with n = 3/2 has been studied by Bilic et al. (2001) in rela- 
tion to the formation of white dwarf stars by gravitational 
collapse. They showed that the quantum pressure can reg- 
ularize the dynamics at small-scales. 
For an isothermal equation of state 



(12) 



(13) 



(14) 



p = p — , 

m 

the effective potential is 

Hp) = — hip, 

m 

and the GP equation reads 
dib h 2 

ih-jfi = -T^ A ^ + + 2k BTeff In 

Interestingly, we note that a nonlinear Schrodinger equa- 
tion with a logarithmic potential similar to Eq. (|14[) has 
been introduced long ago by Bialynicki & Myciclski (1976) 
as a possible generalization of the Schrodinger equation in 
quantum mechanics. 

Remark: for a BEC at T = 0, the pressure arising in 
the Euler equation ([5]) has a meaning different from the 
kinetic pressure of a normal fluid at finite temperature. It 
is due to the self-interaction of the particles encapsulated 
in the effective potential h(p) and not to thermal motion 
(since T = 0). As a result, the pressure can be negative 
(!) contrary to a kinetic pressure. This is the case, in par- 
ticular, for a BEC described by the equation of state (JHJ) 
when the scattering length a s is negative. In terrestrial BEC 
experiments, some atoms like 7 Li have a negative scatter- 
ing length (Fedichev et al. 1996). On the other hand, the 
constant T e ff appearing in Eq. (fT2"j) is just an "effective" 
temperature since it arises from a particular form of self- 
interaction and has nothing to do with the kinetic temper- 
ature (which here is T = 0). In particular, this effective 
temperature can be negative. In that respect, we note that 
linear equations of state p = ape 2 with a < have been 
introduced heuristically to account for the accelerated ex- 
pansion of the universe (see the review of Peebles & Ratra 
2003 and Sec.©. 

3. Newtonian cosmology 

In this first part of the paper, we use the Newtonian cosmol- 
ogy introduced by Milne and McCrea & Milne in 1934. The 
great advantage of the Newtonian treatment is its simplic- 
hjfl Furthermore, the equations are identical with those 
derived using the theory of general relativity, provided the 
pressure is negligible in comparison with the energy density 
pc 2 where c is the speed of light (we shall go beyond this 
limitation in the second part of the paper). This makes this 
approach attractive. Furthermore, its simplicity allows us 



3 It is surprising to realize that Newtonian cosmology was de- 
veloped after, and was influenced by, the cosmological models 
based on Einstein's theory of general relativity. It could have 
been developed much earlier. As Milne writes: "It seems to have 
escaped previous notice that whereas the theory of the expand- 
ing universe is generally held to be one of the fruits of the theory 
of relativity, actually all the phenomena observable at present 
could have been predicted by the founders of mathematical hy- 
drodynamics in the eighteenth century, or even by Newton him- 
self". 



to introduce novel ingredients such as the quantum pres- 
sure which is derived from the classical Gross-Pitaevskii 
equation. The validity and limitation of Newtonian cosmol- 
ogy have been discussed by various authors such as Layzer 
(1954), McCrea (1955), Callan et al. (1965) and Harrison 
(1965). 

We recall the basics of Newtonian cosmology. We con- 
sider a spatially homogeneous solution of Eqs. (@|-([7]) of the 
form 



p(r,t) =Pb(t), u(r,i) 



(15) 



where a(t) is the scale factor and H = a/ a is the Hubble 
"constant" (in fact a function of time) . The Euler equations 
reduce to 



dp b a 

-TT +3p b - 

at a 



-r = -V* 6 . 
a 



0. 



(16) 



(17) 



The pressure p and the quantum potential Q do not appear 
in the theory of the homogeneous model since they enter 
Eq. (|5|) only through their gradients. Equation (fTB"]) leads 
to the relation 



p b a 



1. 



(18) 



which corresponds to the conservation of mass. Taking the 
divergence of Eq. (fT7f and using the Poisson equation d6]), 
we obtain the cosmological equation 



d 2 a 4 _ A 
— = - r G Pb a+-a. 

Using Eq. (fTB")) . its first integral is 
1 



da V 
dt) 



iirGpb + A) a" — k, 



(19) 



(20) 



where K is a constant of integration. Equations (fT5|) - ([2T)| 
are the Newtonian equations of an isotropic and homoge- 
neous universe. They coincide with the equations derived 
by Friedmann (1922,1924) for A = and k = ±1 and by 
Einstein & de Sitter (1932) for A = and k — from the 
theory of general relativity when the pressure is small com- 
pared to the energy density pc 2 . In that case, k is the cur- 
vature constant and space is flat (k = 0), elliptical (k = 1) 
or hyperbolic (k = — 1)Q. 

The Einstein (1917) static universe corresponds to 



0=1, 



Pb = 



A 



k = A. 



(21) 



However, this universe is unstable against perturbations in 
a (Eddington 1930, Harrison 1967). The Einstein-de Sitter 
(EdS) universe corresponds to A = and k = 0. In that 
case, Eqs. (fl"§|) and (|2T))) reduce to 



-7rGp b a, a 2 = ^7rGp b a 2 . 



This yields 
oat" 3 . 



a 6t 



Pb 



6i:Gt 2 ' 



(22) 



(23) 



4 Note that the constant k in Eq. (|2Up can be set to unity by 
a suitable normalization of the parameters. 
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Note that both inflationary theory (Guth 1981) and obser- 
vations favor a flat universe (k = 0). In that case, Eq. (j20|) 
can be written H 2 = ^TrG(p b + pa), where p\ = A/8irG 
is the dark energy density, and it gives a relationship 
between Hubble's constant and the total density of the 
universe. With the present-day value of the Hubble con- 
stant H = 2.273 10~ 18 s _1 , the Einstein-de Sitter model 
(k = A = 0, p <C pc 2 ) leads to an age of the universe 
t = 2/3H Q ~ 9.3 billion years while the ACDM model, 
which is the standard model of cosmology, predicts 13.75 
billion years. The effect of dark energy will be considered 
in Sees. El and |U 

Remark: We can obtain equations (|18l) - ([2"0)) directly if 
wc use the "naive" picture that the universe is a uni- 
form sphere of radius a(t), density pb{t) and mass M. 
The conservation of mass M = ^wp b a 3 leads to Eq. (|18p. 
On the other hand, applying the Gauss theorem at the 
border of the sphere, Newton's equation can be written 
d 2 a/dt 2 = -GM /a 2 + Aa/3 leading to Eq. (HHJ). This is the 
equation of motion of a Active particle of position a in a po- 
tential V(a) = —GM/a — Aa 2 /6. Its first integral is given 
by Eq. (f2T)]) where E = —k/2 can be regarded as the energy 
of the Active particle. The case E = (EdS universe) cor- 
responds to the situation where the velocity of the Active 
particle is equal to the escape velocity. In this naive picture, 
a may be interpreted as the "radius" of the universe. In fact, 
a better derivation (which does not assume a finite sphere) 
is to apply Newton's equation to an arbitrary fluid particle 
located in r and write d 2 r/dt 2 = -G(4:Trp b r 3 /3)/r 2 + Ar/3. 
Setting r = a(f)x and dividing by x yields Eq. (|T9"|) . 



4. Quantum Euler-Poisson system in an expanding 
universe 

We shall now study the instability of the homogeneous 
background and the growth of perturbations that ulti- 
mately lead to the large-scale structures of the universe. 
We shall work in the comoving frame (Peebles 1980). To 
that purpose, we set 



a(i)x, 



a 

u = — r 

a 



(24) 



where v is the peculiar velocity. For the moment, we al- 
low arbitrary deviations from the background flow. Let us 
first write the Poisson equation ([5]) in the comoving frame. 
Integrating Eq. ([T7|) and using Eq. ([T9")) . wc find that the 
background gravitational potential is 



1 a 2 1 •■ 2 2 n i ,\ 2 A 2 
---r =--aax = -irGp b {t)r - -r . 
I a I 3 6 



(25) 



This result can also be obtained by integrating the Poisson 
equation for a homogeneous system or by using the Gauss 
theorem. If we introduce the new potential <j) = <f> — <!>;,, i.e. 



1 2 
(/> = $ + — aax , 

the Poisson equation ((6]) becomes 
A(t> = ATrGa 2 {p-p b ) 7 



(27) 



where the Laplacian is taken with respect to x. Now, fol- 
lowing Peebles (1980) and taking the quantum pressure into 



account, we find that the hydrodynamic equations (U) and 
([5]) can be written in the comoving frame as 



dp 3d 1 . . 

-£ + — P+-V- pv =0, 

at a a 



(9v 1 . . d 1 1 1 

— + -(v • V)v + -v = Vp - -V<f> 

at a a pa a 



(28) 



(29) 



2m 2 a 3 \ y/p J 

It is convenient to write the density in the form 

p = Pb(t)[l + 5(x,t)], (30) 

where p b oc 1/a 3 and <5(x, t) is the density contrast (Peebles 
1980). Substituting Eq. §5Q into Eqs. (J2Z1)- (EHJ) , we ob- 
tain the quantum barotropic Euler-Poisson system in an 
expanding universe 



§ + iv.((l + *)v) = 0, 



dv 1 , a 1 _ 1 , 

— + -(v • V)v + -v = Vp - -V(/> 

at a a pa a 



(31) 



/ AVI 



2m 2 a 3 V VTT5 



Ac/) = 4:TrGp b a 2 S. 



(32) 



(33) 



For h = p = 0, we recover the usual Euler equations of a 
cold gas (T = 0) in an expanding universe (Peebles 19800. 
For a barotropic equation of state p = p(p), the Euler equa- 
tion (|3"2"]l can be rewritten 



dv 1 . _. d 

+ _ v V v+-v = - 7 - 
ot a a (1 



c 2 1 
+- o)a a 



/AVI 



2m 2 a 3 V VTTS 



(34) 



where c 2 = p'(p) = p'(p b (l + 5)) is the square of the veloc- 
ity of sound in the evolving system. It generically depends 
on position and time. For an isothermal equation of state 
c s = hsT/m is a constant and for a polytropic equation 
of state c 2 = K-fp~ f ~ 1 . For a standard BEC with a quartic 
self-interaction, described by the equation of state ©, the 
square of the velocity of sound is 



4ira s h, p 



and we obtain the system of equations 



85_ I 

(26) dt + a 



V • ((1 + £)v) = 0, 



(35) 



(36) 



If the universe is a classical cohisionless fluid, the evolution 
of this fluid is fundamentally described by the Vlasov equation 
(Gilbert 1966). In that case, the hydrodynamic equations with- 
out pressure (h — p = 0) are based on some approximations 
(Peebles 1980) and they are not valid for all times. Alternatively, 
if dark matter is a BEC, the hydrodynamic equations (|31|) - (|33p 
are rigorously equivalent to the GPP system for all times. 
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<9v 1 

dt a 



(v • V)v - 
1, 



4ira s h pb 



VS 



— V0- 
a 



2m 2 a 3 \ y/T+~S 



A(j) = 4irGp b a 2 5. 



5. Linearized equations 



(37) 
(38) 



Of course, 5 = <fi = and v = is a solution of the quan- 
tum Euler-Poisson system ([3~Tj) - (|3"3"|) in the comoving frame, 
corresponding to the pure background flow (|15[) . However, 
this solution may be unstable to some perturbations. If we 
consider small perturbations i5 < 1, ^ < 1, |v| <C 1 and 
linearize the foregoing equations we obtain 



85 1 
at a 



0, 



(39) 



<9v a 

757 + _v 

at a 



--c 2 V6- 
a a 



Am 2 a 3 



V(A<5), 



A(f> = AnGp b 8a 2 



(40) 
(41) 

where c 2 s = p'(pb(t)) now denotes the square of the velocity 
of sound in the homogeneous background flow. It is just 
a function of time. Taking the time derivative of Eq. (|3"9"|) 
multiplied by a, the divergence of Eq. (|4T)|) . and using Eq. 
(|4"Tj) . these equations can be combined into a single equation 
governing the evolution of the density contrast 



d 2 S 



a 85 



dt 2 + 2 adt 



= -^A5 + 4nGp b S 



tr 



-A 2 5. 



For h = 0, we recover the equation first derived by Bonnor 
(1957). In his case, the pressure p is a kinetic pressure. 
Expanding the solution in Fourier modes of the form 
<5(x, t) = <!>k(t)e lk ' x , we obtain 



a \ 4m z a^ 



2 1.2 



cjk 



4irGp b )S = 0, 



(43) 



where, for brevity, we have noted 5(t) for 5k(t). For a stan- 
dard BEC, using c 2 = A-Ka s fi 2 pb/m 3 , the foregoing equation 
can be rewritten 



a \ 4m z a^ 



4Tra s h pbk 2 



iirGpb )S = 



(44) 



6. Quantum Jeans length 

6.1. Expanding universe 

In the non-interacting case a s = c s = 0, the equation for 
the density contrast reduces to 



a- ( h 2 k A 
5 + 2-5+1 ——r 



- inGpb (5 = 



(45) 



From this relation, we can define a time-dependent quan- 
tum Jeans wavenumber 



/ IGirGpb 



2„4\ 1/4 



(46) 



Note that the proper quantum Jeans wavenumber ob- 
tained by writing 6 cx e lk * r is k* Q = (16irGp b m 2 /h 2 ) 1 / i . 

Recalling Eq. (fT8)l . we can write kQ = kqo 1 / 4 where 
kq = (16nGpba 3 m 2 /h 2 ) 1 ^ is a constant. The quantum 
Jeans length Xq = 2n/kQ decreases with time like a -1 / 4 so 
that, in the comoving frame, the system becomes unstable 
at smaller and smaller scales as the universe expands (note 
that, on the contrary, the proper quantum Jeans scale in- 
creases with time like a 3 / 4 ). In the cold, non-quantum uni- 
verse, there is no Jeans length: all the scales are unstable. 
Therefore, quantum effects can stabilize the system at small 
scales and avoid the density cusps (Hu et al. 2000). 

In the TF limit where the quantum potential can be 
neglected, the equation for the density contrast reduces to 



5 + 2-5 

a 



- 4nGp b )5 = 



(47) 



From this relation, we can define a time-dependent classical 
Jeans wavenumber 



/ A-KGpbO 2 



1/2 



(48) 



The proper Jeans wavenumber is kj = (AnGpb/cl) 1 ^ 2 . The 
evolution of kj(t) depends on the equation of state. For 
an isothermal equation of state for which c 2 = ksT/m, 
recalling Eq. (|18[) . we can write kj = KjaT 1 ! 2 where 
kj = (A-KGpbcfim/kBT) 1 / 2 is a constant. The classical 
Jeans length Aj = 2n/kj increases with time like a 1 / 2 so 
that, in the comoving frame, the system becomes unstable 
at larger and larger scales as the universe expands. For a 
(42) polytropic equation of state for which c 2 = KjpJ -1 , we can 



write k 



J 



Kja 



(3 7 - 



4 )/ 2 where kj = [4irG(p 



6 a 



3^2-7 



/Kj] 



1/2 



is a constant. The classical Jeans wavelength behaves like 
. For 7 < 4/3, the Jeans wavelength grows 



A 



(4-3 7 )/2 



j (x a 

with time and for 7 > 4/3, it decreases with time. For 
7 = 4/3, the Jeans wavenumber is constant in time, i.e. 
kj = kj. Finally, for a standard BEC, c 2 = 4ira s h pb/m 3 , 
and the Jeans wavenumber can be written 



k.,= 



Gm 3 a 2 
a s h 2 



1/2 



(49) 



It is independent on the density and can be written kj = 
Kja with kj = (Gm 3 / ash 2 ) 1 / 2 . The Jeans length decreases 
like A j cx 1/a. 

Remark: the study of the Jeans instability in an ex- 
panding universe exhibits a critical value of the polytropic 
index 'y cr it =4/3 (i.e. n cr i t = 3). It is interesting to note 
that the same critical index arises when one studies the 
dynamical stability of spatially inhomogeneous polytropic 
spheres with respect to the barotropic Eulcr-Poisson sys- 
tem. It has been established that a polytropic star is stable 
for 7 > 4/3 and unstable otherwise (Binney & Tremaine, 
1987). Surprisingly, the same index appears in the stability 
analysis of a spatially homogeneous polytropic gas in an ex- 
panding universe. Note that the static study of Jeans (see 
the following section) docs not directly exhibit a critical 
value of the polytropic indc>0. 



The critical value jcrit = 4/3 only appears when one intro- 
duces the Jeans mass Mj ~ pXj tx pt 37-4 '/ 2 , b u t n does not 
play any particular role in the stability analysis of the homoge- 
neous fluid. 
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6.2. Static universe 

If we assume that Eq. (1431) remains valid in a static universe, 
and take a = 1, we get 



5 + 



( 



2 1.4 



ffk 



V 4m 2 



+ ctk 



2 ;„2 



47rGp 6 ) (5 = 0. 



(50) 



Writing the perturbation in the form S(t) oc e we ob- 
tain the dispersion relation 



ft 2 fc 4 
Am 2 



c 2 s k 2 - 4irGp b 



(51) 



This dispersion relation has been studied in our previous 
paper (Chavanis 2011). The generalized Jeans wavenumber 
is given by 



k 2 = 



2m 2 



'c* + 



AirGh 2 



Pb 



(52) 



In the non- interacting case (a s = c s — 0), we recover 
the quantum Jeans wavenumber (|46[) with a = 1 and in 
the TF limit, we recover the classical Jeans wavenumber 
l|48p with a=l. The dispersion relation can be written 
u! 2 /AirGpb = fc 4 /kq + k 2 /k 2 — 1 and the generalized Jeans 
wavenumber fc 2 = (fc 4 /2fc 2 )[±(l + Akj/k Q )^ 2 - 1] with 
+ when c 2 > and — when c 2 < 0. For A < A c , w is 
real and the perturbation oscillates with a pulsation u>; 
for A > A c , u> is purely imaginary and the perturbation 
grows exponentially rapidly with a growth rate 7 = v 1 —u) 2 . 
Jeans' stability criterion is therefore A < A c . For c 2 > 0, 
the maximum growth rate corresponds to fc = (infinite 
wavelengths) and is given by 7 = yJAirGpb- On the other 
hand, considering a BEC with negative scattering length, 
it is found (Chavanis 2011) that the maximum growth rate 
corresponds to fc* = (& j n\a s \pb/m) 1 ' 2 and is given by 



7* 



'l67T 2 a 2 h 2 p 2 



AirGpb- 



(53) 



Note that fc* = (fc 4 ? /2|fc 2 T |) 1 / 2 and 7* = y/HrUp(l + 
fc g/ 4fc j) 1/2 - For > (Gm 4 /4irh 2 p b )V 2 , we find that 
7* ~ 4ir\a s \frpb/m 2 3> y/AwGpb- Therefore, an attractive 
short-range interaction (a s < 0) increases the growth rate 
of the Jeans instability. 

Jeans' classical analysis (Jeans 1902, 1929) suffers from 
the defect that in general there is no initial stationary state 
that is in a uniform non-rotating fluid. Therefore, using 
Eq. (|50[) in a static universe is referred to as the "Jeans 
swindle"Q. As noted by Bonnor (1957), the Jeans procedure 



7 Kiessling (2003) provides a vindication of the "Jeans swin- 
dle". He argues that, when considering an infinite and homo- 
geneous distribution of matter, the Poisson equation must be 
modified so as to correctly define the gravitational force. He 
proposes to use a regularization of the form A$ — &□<!> = AirGp 
where ko is an inverse screening length that ultimately tends to 
zero (ko -> 0), or a regularization of the form A 4? = 4nG(p — ~p) 
where p is the mean density. In his point of view, this is not a 
swindle but just the right way to make the problem mathemat- 
ically rigorous and have a well-defined gravitational force. 



is only valid in the static Einstein universal- However, in 
this last case, the background density of the universe is 
Pb = A/AirG. The justification does not apply to a uniform 
mass of gas of different density. Furthermore, the Einstein 
universe is strongly unstable. Therefore, this justification is 
not valid and it is necessary to develop the Jeans instability 
analysis in an expanding universe. 

Remark: Jeans' analysis could be valid if the growth 
rate of the instability were much larger than the rate of the 
expansion of the universe. According to Eq. (|2"2")l. the rate 
of expansion of the universe is H = a/a = (STtGp/S) 1 ' 2 . On 
the other hand, when c 2 > 0, the maximum growth rate of 
the Jeans instability is 7 = \fAnGp. They are exactly of 
the same order which leads to the usual conclusion that 
Jeans' instability analysis is not valid (see, e.g., Weinberg 
1972). However, when c 2 < 0, the maximum growth rate, 
given by Eq. ([53)1 can be much larger. In that case, the 
expansion of the universe could be ignored and the static 
Jeans instability analysis could be valid. 



7. Solution in the Einstein-de Sitter universe 

The evolution of the perturbations is more complicated to 
analyze in an expanding universe than in a static universe 
(and turns out to be very different). For simplicity, we con- 
sider the case of an Einstcin-de Sitter universe. In that case, 
it is possible to obtain analytical solutions of the linearized 
equation (|43[) for the density contrast in some particular 
cases. Measuring the evolution in terms of a instead of t 
and using Eq. (|2"2"j) . Eq. (|4"3"]l is transformed into 



d 2 d 3_dS_ 

da 2 2a da 



2;,4 



Tik 



2a 2 \16irGp b m 2 a 4 AnGp b a 2 



-1)5 = 0. 



(54) 



For h — c s = 0, we recover the case of a cold classical 
universe in which the pressure is zero and quantum effects 
are neglected. Equation (j54|) reduces to 



d 2 5 + 3_dS_ _ 3 5 
da 2 2a da 2a? 

and its solutions are 



0. 



5_ oc a- 3 ' 2 . 



(55) 



(56) 



We note that the growth of the perturbations is algebraic in 
an expanding universe while the "naive" Jeans instability 
analysis in a static universe predicts an exponential growth. 
Therefore, the results of the two analysis are very different. 
As mentioned long ago by Bonnor (1957) and others, this 
(slow) algebraic growth may be a problem to form large- 
scale structures sufficiently rapidly. 

Assuming a polytropic equation of state p = Kp 1 for 
which c 2 = Kjp2 1 and introducing the variables of Sec. 



8 As is well-known, Einstein (1917) introduced a cosmological 
constant in the equations of general relativity in order to recover 
a homogeneous static universe. He also considered, as a preamble 
of his paper, a modification of the classical Poisson equation in 
the form A$ — ko$> = AivGp because he (incorrectly) believed 
that, in the Newtonian world model, the cosmological constant 
is equivalent to a screening length (see Spiegel 1998, Kiessling 
2003 and Chavanis & Delfini 2010 for historical details). 
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16.11 we can rewrite Eq. in the form 




d 2 8 3_dd_ 

da 2 2a da 



Making the change of variables 8(a) 
the differential equation for / is 



lU = 0. 



(57) 



f(a)/a, we find that 



da 2 



2a da 



1 



3fc 4 



3fc 2 



2k 2 o 37 



- - 1 / = 0. (58) 



We have not been able to find the general solution of this 
equation in terms of simple functions so that we shall con- 
sider particular cases. Some of the solutions presented be- 
low have already been discussed in the literature (see, e.g., 
Harrison 1967) but we shall give more detail and put em- 
phasis on the solutions describing self-gravitating BECs 
that were not discussed previously. 

7.1. The TF approximation 

In the TF approximation where the quantum potential can 
be neglected, Eq. (|58j) reduces to 



da 2 



2a da 



1 



3fc 2 



2K 2 a^- A 



- 1 / = 0. 



(59) 



This is a particular case of the equation studied by Savcdoff 
& Vila (1962) in relation to the work of Bonnor (1957)0. It 
can be solved in terms of Bessel functions (see Appcndix[2]) • 
For 7 ^ 4/3, f(a) is given by Eq. (|A.7[) and the evolution 
of the density contrast is 



5(a) oc 



aJZ J ±z 



y/6 k 4-3 7 

a 



37 Kj 



(60) 



Using the asymptotic expansions of the Bessel functions, 
we note that the density contrast is oscillating for k 
kj(a) and growing for k <C kj(a) (on the timescale over 
which these inequalities are fulfilled). Therefore, wave- 
lengths which are long with respect to \j(a) behave like 
in a cold universe whereas for the short wavelengths the 
perturbations have an oscillatory behavior. In this sense, 
the results are similar to those obtained with the classi- 
cal Jeans analysis. They are in fact different because (i) 
the Jeans length varies with time, (ii) the oscillating so- 
lutions can be growing or decaying and (iii) the growth is 
algebraic instead of exponential. In fact, the results more 
crucially depend on the value of the polytropic index 7 than 
on the Jeans length Aj. We must distinguish four cases, (i) 
7 < 4/3: for a — > the perturbations grow like in a cold gas 
~ a) and, for a — > +00, they undergo damped oscilla- 
tions (they decay like I/a* 5 " 37 )/ 4 ); (ii) 4/3 < 7 < 5/3: for 
a — >• the perturbations diverge like l/o/ 5-37 ")/ 4 while os- 
cillating implying that we must start the stability analysis 
at at > 0. For a — > a;, the perturbations undergo damped 
oscillations (they decay like l/a^ 5-37 ^ 4 ) and, for a — > +00, 
they grow like in a cold gas (S + ~ a); (iii) 7 = 5/3: for 
a — > the perturbations oscillate and, for a — > +00, they 
grow like in a cold gas (8+ ~ a); (iv) 7 > 5/3: for a — > 



In this section, we allow p to have the interpretation of an 
ordinary pressure, i.e. we consider not only a gas of BECs but 
also a more conventional gas with a kinetic pressure. 



the perturbations undergo growing oscillations (they grow 
like a/ 37-5 -*/ 4 ) and for a +00 they grow like in a cold gas 
(8+ ~ a). Considering the regime a <C 1, the perturbations 
grow for 7 < 4/3, decay for 4/3 < 7 < 5/3 and grow for 
7 > 5/3. Considering the regime a 3> 1, the perturbations 
decay for 7 < 4/3 and grow for 7 > 4/3. For 7 < 4/3, the 
perturbations start to grow but finally decay. In a sense, 
the system is asymptotically stable. However, the growth 
of the perturbations in the initial stage can trigger nonlin- 
ear effects that may induce instabilities. 

For 7 = 4/3 (n = 3), the density contrast behaves like 



8(a) oc a 



(61) 



This polytropic index corresponds to a gas of photons at 
temperature T or to a gas of relativistic fermions at zero 
temperature. In that case, the Jeans length is independent 
on time: kj = kj. For k > (25/24) 1 / 2 fcj the density con- 
trast behaves like 



8(a) oc 



1 



cos 




(62) 



It diverges for a — > implying that we must start the sta- 
bility analysis at ai > 0. For a > a%, the perturbations 
decay like a -1 / 4 while making oscillations. For kj < k < 
(25/24) 1 / 2 fcy, the perturbations decay algebraically with- 
out oscillating. For k < kj, the perturbations grow alge- 
braically (more precisely 8+ grows and 8- decays). For 
k — > 0, the perturbations behave like in a cold gas (5+ 
grows like a and <5_ decays like a -2 / 3 ). This situation is 
relatively close to the classical Jeans stability analysis but 
(i) the growth of the perturbation is algebraic instead of 
exponential; (ii) the oscillatory solutions decay with time; 
(iii) there is a small interval kj < k < (25/24) 1//2 fcj that 
has no counterpart in the static Jeans analysis. 

Let us consider particular equations of state with a 
physical meaning. The index 7 = 1 corresponds to an 
isothermal universe. The density contrast is given by 



8(a) oc 



1 



,1/4 



J±5 ( V6 



,1/2 



For the index 7 = 5/3 (n = 3/2), we obtain 



8(a) oc 



1 



,1/4 ' J ±i 



V6- 



1 



Kj a 



1/2 



(63) 



(64) 



This polytropic index corresponds to a gas of non- 
relativistic fermions at zero temperature. For the polytropic 
index 7 = 2 (n = 1), we have 



8(a) oc 




(65) 



This polytropic index corresponds to a BEC with quartic 
self- interaction in the TF approximation (see Sec. [5]). In 
that case, kj = (Gm 3 /ash 2 ) 1 / 2 . The previous expression 
assumes that a s > 0. In the case where a s < 0, we obtain 



8(a) oc 




(66) 



<s 



P.H. Chavanis: Growth of perturbations in an expanding universe with Bose-Einstein condensate dark matter 



with kj = (Gm 3 /\a s \ri 2 )^ 2 . For a - 
tions behave like in a cold gas. For o 

5(a) oc a 1 / 4 e v ^~~ — s- +00. 



+00, the perturba- 
0, we find that 

(67) 



In a sense, this exponential divergence demonstrates that 
an attractive self-interaction accelerates the growth of the 
perturbations. In that case, we can neglect the expansion 
of the universe and we are led back to the static study of 
Sec.El 

Some curves representing the evolution of the density 
contrast 5(a) in the different cases listed above are repre- 
sented in Figs. Q][S] for illustration. 




Fig. 1. Evolution of the perturbation 5(a) in an isothermal 
universe (7 = 1). This corresponds to the case 7 < 4/3. We 
have taken k — kj. 



Y=4/3 



1 < k/kj < (25/24) 




k/kj > (25/24) 



0.5 



Fig. 2. Evolution of the perturbation 5(a) in a photonic 
or in a relativistic fermionic universe (7 = 4/3). We have 
represented the solution <5 + . 



7.2. The non-interacting BEC case 

In the non-interacting case a s = c s = 0, which corresponds 
to a BEC universe without self-interaction, Eq. (|58[) reduces 
to 



da' 2 2a da 



1 



3fc 4 



- 1 / = 0. 



(68) 



0.5 



-0.5 




y= 1.55 



2 3 4 

a 



Fig. 3. Evolution of the perturbation 5(a) in a polytropic 
universe with 7 = 1.55. This corresponds to the case 4/3 < 
7 < 5/3. We have taken k = kj. 




Fig. 4. Evolution of the perturbation 5(a) in a non- 
relativistic fermionic universe or in a BEC universe without 
self- interaction (7 = 5/3). We have taken k = Kj. 




Fig. 5. Evolution of the perturbation 5(a) in a BEC uni- 
verse with quartic self-interaction in the TF limit (7 = 2). 
This corresponds to the case 7 > 5/3. We have taken 
k = kj. 



Comparing this equation with Eq. (|59[). we see that the 
dependance in a is the same as for a polytrope 7 = 5/3. 
Therefore, f(a) is given by Eq. (|A.7|) with A = 3k 4 /2n^ 
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and a = — 1/2 yielding 



The equation for /(a) is 



5(a) oc 



V6 



k' 2 



K n a 



1/2 



(69) 



For a — > the perturbations oscillate and for a — > +oo they 
behave like in a cold gas (5+ grows like a). In a sense, a cold 
bosonic universe behaves similarly to a cold non-rclativistic 
fcrmionic universe (compare Eqs. (|69[) and ([64])). However, 
the dependence in k is different (k instead of k) and the 
expression of the Jeans length also is different (compare 
Eqs. g6l) and ((48]) ). 

7.3. The non-relativistic fermionic case 7 = 5/3 

A non-relativistic gas of fcrmions at T = is usually 
described by a polytropic equation of state p = Kp 1 
with index 7 = 5/3 (n = 3/2) and polytropic constant 
K = (1/5)(3/8tt) 2 / 3 /i 2 /to 8/3 (Chandrasekhar 1939). This 
equation of state arises from the Pauli exclusion principle. 
Such a semi-classical description is valid in the TF approx- 
imation which is exact when TV — > +00. In some cases, it 
can be relevant to go beyond the TF approximation and 
take into account the quantum pressure arising from the 
Hciscnberg principle. This can regularize the dynamics at 
small scales as shown by Bilic et al. (2001) in their study 
of the formation of fcrmion stars by gravitational collapse. 
In our problem, this general situation is described by Eq. 
([57]) with 7 = 5/3. This yields 



da 2 



3 dS 3 
2a da 2a 2 



The equation for f(a) is 

d 2 f 1 df if 3fc 4 3fc 2 



da 2 2a da a 2 \ 2K%a 2n 2 a 



16 = 0. 



1 / = o. 



(70) 



(71) 



As noted in Sec. 17.21 the two terms have the same depen- 
dence on a. Therefore, f(a) is given by Eq. (|A.7[) with 



A- — + — a- -- 

2/"C j 2 

We obtain 

1 T ( [6k~ 4 6k 2 1 \ 

For a — > the perturbations oscillate and, for a — > +00, 
they behave like in a cold gas. 



(72) 



(73) 



7.4. The relativistic fermionic case 7 = 4/3 

A relativistic gas of fcrmions at T = is usually de- 
scribed by a polytropic equation of state p = Kp 1 with 
index 7 = 4/3 (n = 3) and polytropic constant K = 
(l/4)(3/87r) 1 / 3 /ic/?n 4 / 3 (Chandrasekhar 1939). If we go be- 
yond the TF approximation and take the quantum pres- 
sure into account, the evolution of the density contrast is 
described by Eq. <j57j) with 7 = 4/3. This yields 



d 2 S 3_dS_ 

da 2 2a da 




- 1 \S = 0. 



(74) 



da 2 2a da a 2 \ 2kqO ' 2k 2 



/ = 0. 



(75) 



Its solution is given by Eq. (|A.13[) and the evolution of the 
density contrast is 



5(a) OC : ; , J 

v ' a 1 / 4 ±A 



25 — 24- 



t h/6 



k 2 1 

k 2 q a 1 / 2 



(76) 



We note that the classical Jeans scale appears in the in- 
dex of the Bessel function while the quantum Jeans scale 
appears in the argument of the Bessel function. 

7.5. The general case (asymptotics) 

We now treat the general case, taking into account quan- 
tum pressure and polytropic pressure (we assume a s > 0). 
Since we cannot solve Eq. (|57p analytically, we just give 
asymptotic results. 

Let us first consider the limit a — >• 0. (i) If 7 < 5/3, 
the polytropic pressure is negligible in front of the quantum 
pressure and the perturbations oscillate (see Sec. I7.2|) . (ii) If 
7 = 5/3, the polytropic pressure and the quantum pressure 
are of the same order, and the perturbations oscillate (see 
Sec. l7.3[) . (hi) If 7 > 5/3, the quantum pressure is negligible 
in front of the polytropic pressure and the perturbations 
undergo growing oscillations 5 oc a^ 37-5 )/ 4 (see Sec. 17. ip . 
Therefore, the perturbations oscillate for 7 < 5/3 and grow 
for 7 > 5/3. The quantum pressure has a stabilizing role 
for 7 < 4/3 (see Sec. EH . 

Let us now consider the limit a — > +00. (i) If 7 < 4/3, 
the quantum pressure is negligible in front of the poly- 
tropic pressure and the perturbations undergo damped os- 
cillations and decay like l/a^" 37 )/ 4 (see Sec. O]). (ii) If 
7 = 4/3, the quantum pressure is negligible in front of 
the polytropic pressure and we are led back to the critical 
case of Sec. 17.11 for k > (25/24) 1 / 2 fcj the perturbations 
decay like a -1 / 4 while making oscillations, for kj < k < 
(25/24) 1 / 2 fcj the perturbations decay algebraically with- 
out oscillating and for k < kj the perturbations grow alge- 
braically, (iii) For 7 > 4/3, the perturbations behave like 
in a cold gas and grows like 5 + ~ a — > +00. Therefore, the 
perturbations decay for 7 < 4/3 and grow for 7 > 4/3. 

8. Bose-Einstein condensate universe 

8.1. Newtonian cosmology with pressure 

We now take into account pressure effects (coming from 
special relativity) in the evolution of the cosmic fluid. We 
base our study on the set of hydrodynamic equations 

^+V-(pu) + 4v-u = 0, 



(77) 
(78) 

(79) 



where c is the velocity of light. For c — » +00, we recover the 
classical Euler-Poisson system (Binney & Tremaine 1987). 



— + u-V u = -V" V$ ' 

ot p + 

A$ = 4^G (p+%\ - A, 
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These equations were introduced by McCrea (1951). His 
initial continuity equation, which contained an incorrect 
pressure gradient term, was corrected by Lima et al. (1997) 
and wc took this correction into account in writing Eq. (1771) . 
These equations lead to the correct relativistic equations for 
the cosmic evolution. Indeed, assuming a homogeneous and 
isotropic solution of the form p(r,t) = pb(t), p(v,t) = pb(t) 
and u(r,i) = (d/a)r, Eqs. ([77 |) -([79 l reduce to 



^fl + 3" ( + El 

dt a V c 2 



= 0, 



a 
a 



4ttG 



Pb 



3pb 

„2 



A 

3' 



8nG 



-Pb 



K 



A 
I' 



(80) 



(81) 



(82) 



These are precisely the Fricdmann equations that can be de- 
rived from the theory of general relativity (Weinberg 1972). 
The perturbation theory based on Eqs. (| 7 7 1) - (179]) has been 
discussed by Reis (2003) who showed that, depending on 
the equation of state, the results may agree or disagree 
with the general relativistic approach. We shall leave this 
problem open since it is not our purpose here to investigate 
the validity of these equations in detail. Furthermore, our 
main results will concern the evolution of the cosmic fluid 
governed by the Friedmann equations. Since the Friedmann 
equations can be derived from the theory of general relativ- 
ity, the validity of these equations is not questioned. Note 
that equations different from Eqs. ([TT|) - (|79|) have been con- 
sidered by Pace et al. (2010) and Harko (2011). Their equa- 
tions involve an additional term pu/c 2 in the Eulcr equation 
(|78p . However, this term apparently leads to equations for 
the cosmic evolution that are different from the Friedmann 
equations, so that this term will not be considered here. 

8.2. The cosmic evolution of a BEC universe 

We shall assume that the dark matter is a BEC with quartic 
self-interaction^ described by the barotropic equation of 
state JH]). For convenience, we write this equation in the 
form 



2 „2 



p = kc p 



(83) 



where k = 2Tra s h /m c is a constant that can be positive 
(a s > 0) or negative (a s < 0). In a first step, we concen- 
trate on the dark matter component and neglect radiation, 
baryonic matter and dark energy (in particular we take 
A = 0). These additional terms will be considered later. 
For the equation of state ([8"3"]) . the Friedmann equation (|5D"|) 
becomes 



d Pb , o h ft \ 
— + 3-p b {l + kpb) 
dt a 



0. 



This equation can be integrated into 
A 

Pb = ■ 



kA 1 



(84) 



(85) 



where A is a constant. If we make the physical requirement 
that equation (|85[) has solutions for large a, then we must 



In this part of the paper, we neglect the quantum pressure. 



impose A > 0. In that case, pb ~ A/a 3 for a — > +00 which 
returns Eq. (fT8j) . When k > 0, which is the case previously 
considered by Harko (2011), the density exists only for a > 

= {kA) 1 ' 3 . For a — > a*, pb — > +00. When k < 0, the 
density is defined for all a and has a finite value pb = 1/ 1 k\ 
when a — > 0. 

Combining Eqs. (|8Tj) and ((83)) . we obtain 



a 
a 



AnGpb 



(l + 3kpb). 



(86) 



When k > 0, the universe is always decelerating (d < 0). 
When k < 0, the universe is accelerating (d > 0) for pb > 
p c = l/3|fc| and decelerating (d < 0) for pb < p c - Using 
Eq. (|85|) . this can be re-expressed in terms of the radiu^l: 
the universe is accelerating for a < a c = (2A|fc|) 1 / 3 and 
decelerating for a > a c . 

In order to determine the temporal evolution of a(t), we 
shall assume k = (flat space) like in the Einstein-de Sitter 
universe. Combining Eqs. (|82|) and (j85|) . we get 



/8ttGA 



1/2 



kA 



(87) 



We note that for a — > +00, this equation reduces to Eq. 
(f2"2"|) with Eq. (|18p and we recover the Einstein-de Sitter 
solution ((23]) . It is convenient to define a* = (IfcjA) 1 / 3 and 
introduce R = a/a*. In that case, the density is given by 



1 



Pb 



k\ R 3 T 1' 
and Eq. ([57} can be rewritten 
KR 



R = 



(89) 



where the upper sign — corresponds to k > and the lower 
sign + corresponds to k < 0. We have denned the constant 



K 



1/2 



4Gm 3 c 2 ^ 1/2 



3\a s \h 



(90) 



Note that it depends on the mass m of the bosons and on 
their scattering length a s . Introducing the dimensionless 
parameter A = a s /X c = a s mc/h (Chavanis 2011) measur- 
ing the strength of the short-range interactions (A c is the 
Compton wavelength of the bosons), we can rewrite Eq. 
(l90l) in the form 



K 



2 / m 
V3 V^p 



1 



p ' 



(91) 



where Mp = (hc/G) 1 ' 2 is the Planck mass, lp = 
(HG/c 3 ) 1 / 2 is the Planck length and tp = Ip/c is the Planck 
time. 

For k > 0, which is the case previously considered by 
Harko (2011), the solution of Eq. ([89]) is 



\JR 3 - 1 - arctan V R 3 - 1 = -^Kt, 



(92) 



We shall sometimes call the scale factor a the "radius" , hav- 
ing in mind the naive picture of a universe behaving like a homo- 
geneous sphere of radius a. This is of course an abuse of language 
in the context of the theory of general relativity. 
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where the constant of integration has been set equal to zero. 
In this model, the universe starts at a finite time t = (see 
Scc. l8.5l for a revision of this statement) with a finite radius 
R(0) = I and an infinite density pb(0) = oo. For t —> 0, 
R ~ 1 + (3/ 4) 1 / 3 (i i f t) 2 / 3 . The universe is expanding, always 
decelerating, and asymptotically approaches the Einstein- 
de Sitter universe R ~ (3Kt/2) 2 / 3 (see Fig. EJ}. 




Fig. 6. Evolution of the scale factor in a BEC universe with 
k > 0. The dashed line corresponds to the Einstein-de Sitter 
universe (k = 0) that is reached asymptotically. 




Fig. 7. Evolution of the scale factor in a BEC universe with 
k < 0. The dashed line corresponds to the Einstein-de Sitter 
universe (k = 0). We note that a BEC universe with k < 
converges more rapidly towards the EdS solution than a 
BEC universe with k > 0. The bullet (locating the inflexion 
point) corresponds to the time t c at which the universe 
starts to decelerate. 



For k < 0, the solution of Eq. ((8S) is 

r—— ( 1 + VW+l\ 3_ 
V i? 3 + 1 — In I ^ = -Kt, 



(93) 



where the constant of integration has been set equal to 
zero. In this model, the universe starts from t — > — oo 
with a vanishing radius R(— oo) = and a finite density 
p 6 (-oo) = l/|fc|. For R< R c = 2 1 / 3 ~ 1.26, the universe is 
accelerating and for R > R c = 2 1 / 3 , it is decelerating. For 



i? — > +oo, we recover the classical Einstein-de Sitter uni- 
verse- R ~ (3i«/2) 2/3 (see Fig. 0). The change of regime 
between the phase of acceleration and the phase of decel- 
eration manifests itself by an inflexion point (R = 0). The 
time at which the universe starts decelerating is 



-Kt c = \/3 - In 



1 + V3 
V2 . 



1.074. 



(94) 



At t = 0, R(0) ~ 0.7601 and p b (0) = 0.6948/|fc|. For t -> 
— oo, i? —t and the asymptotic expansion of Eq. (|93|) yields 



(95) 



For t <C t c , the universe expands exponentially rapidly, 
decelerates after t c and coincides with the Einstein-de Sitter 
universe for t ^> t c . In the accelerating phase, the rate of 
the expansion is given by Eq. ((90)) . 

Important remark: in our mathematical description of 
the evolution of the cosmic fluid, we have extrapolated the 
solutions ((92)) and (|93|) far away in the past. In the model 
with k > 0, the universe emerges at a primordial time t = 
from a big-bang singularity where the density is infinite. In 
the model with k < 0, the universe has always existed (up 
to t — > — oo) and there is no big bang singularity at t = 0. 
This is due to the negative pressure allowing for the pres- 
ence of an inflexion point that changes the concavity of the 
curve a(t). Of course, this extrapolation to t — > — oo is not 
physically justified since we have ignored important effects 
like the radiation which dominates in the early universe. 
We shall see in Sec. 18.41 that the contribution of the radi- 
ation strongly alters the results at early times in the case 
k < 0. Therefore, our description of the solution ((93)) must 
be considered only on a formal mathematical basis. Still it 
shows that one can construct cosmological models (solution 
to the Einstein equations) without primordial singularity. 

8.3. Other representations 

The previous representation shows that there exists a single 
universal curve a{t) in each case a s > 0, a s = (EdS) and 
a s < 0, provided that the units of time and length are 
appropriately chosen. However, these units depend on a s . 
Alternatively, it may be useful to choose units of time and 
length that are independent on a s and plot the curve a(t) 
for different values of a s . This is the representation that has 
been chosen by Harko (2011) for the case k > and that 
we shall discuss and generalize in this section. 

In Eq. I|85p. the constant A can be determined by the 
present-day density pb = po and the corresponding scale 
factor a = uq. Writing po = ^V( a o — we obtain 



A = 



l + kpo' 



(96) 



We note that po < l/\k\ when k < as we have already 
indicated. The relation (|85)) between the density and the 
scale factor can be rewritten 



q 

Pb = T 



1 



where we have introduced the dimensionless parameter 
kpo 



(97) 



<1 



1 + kp 



(98) 
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We note that q is positive for k > and negative for k < 0. 
Furthermore, q = for k = 0, q — > 1 for k — > +oo and q — > 
— oo for k —1/po. Using these results, the differential 
equation ([57)1 giving the evolution of the scale factor with 
the time can be rewritten 



inGq 
3fc 



1/2 



(-) 



(99) 



From Eq. ([82)1 , we can define the present-day critical density 
by (Pc)o = 3Hq/8ttG where i?o is the present-day Hubble 
constant. If k = 0, which seems to be the case for our 
universe, (p c )o represents the total density including bary- 
onic matter, radiation, dark matter and dark energy. It is 
customary to introduce the present-day density parameter 
= Po/(Pc)o- Then, we can write 



3k 



v 7 ! + kp 



(100) 



Since the denominator depends on k, it must be expressed 
in terms of q (which is our control parameter). Using Eq. 
(|98p. we obtain 



kp = 



<1 



so that Eq. (|99p can be rewritten 
a = H ^U~ (1 - qf' 2 



v^ 3 



The density is 
1-9 



Pb = Po 



(101) 



(102) 



(103) 



The curves giving the scale factor a(t) as a function of 
time for different values of q are plotted in Figs. \8\ and 
[HI The interest of this representation is that it shows how 
the BEC models (|105[) and (|106[) approach the Einstcin-dc 
Sitter model (|10T[) as q — > 0. On the other hand, for larger 
values of q (positive or negative), there can be substantial 
differences between a BEC universe and the classical pres- 
surelcss Einstcin-dc Sitter universe. 



k>0 








/ / 

i,i, 



1/2 

Q -H t 

Fig. 8. Scale factor a(t) as a function of time t for a BEC 
universe with k > 0. The units of space and time are nor- 
malized by present-day quantities that do not depend on 
k. The different curves correspond to q = (EdS, dashed), 
q = 0.01, q = 0.1 and q = 0.5. The straight line corresponds 
to a/ciQ — 1. Its intersection with the curve a(t) defines the 
age of the universe. 



We can make the connection with the notations of Sec. 18.21 
by setting 



ao\1 



1/3 



'I 



1/2 



(104) 



In that case, Eq. (|102[) takes the form of Eq. ([59")) and its 
solutions are given by Eqs. ([9"2")l and ([9"3")l . Returning to the 
notations of this section, we find for < q < 1 that 



7iMvO^"^ arctan ' 



(-) 3 

^ an ' 



(105) 



1 1 

k<0 






/// EdS 

,i,i,i 



1 

1/2 

n H t 



which was previously obtained by Harko (2011). On the Fi S- 9 - Scale factor °(*) as a function of time t for a BEC 



other hand, for q < 0, we obtain 



-75^( A /(-) 3 -9-A/Mln 
Vl - q I V a 



3 



universe with k < 0. The units of space and time are nor- 
malized by present-day quantities that do not depend on 
k. The different curves correspond to q = (EdS, dashed), 
q = -0.01, q = -0.1 and q = -0.5. 



-Vn H t. (106) 



To see that, let us consider the asymptotic behaviors of 
Finally, for q = 0, we recover the Einstein-de Sitter model Eqs. (j!05p and (|106p . For i = 0, we find that 



(107) ^ = i?±M i/3, 
a 



(108) 
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with R + = 1 for k > and R. 

t — > +00, we get 



0.7601 for k < 0. For The age of the EdS universe (q = 0) is 




q) 1/3 (H tf/\ 



(109) 



For < q < 1, the scale factor is asymptotically smaller 
than in an EdS universe (corresponding to q = 0). Since 
a BEC universe with k > initially starts with a radius 
a(0) > 0, this implies that the curves must cross each 
other at some point (such a crossing is shown in Fig. [5] 
for q = 0.5). Of course, the crossing point occurs at larger 
and larger times as q — > 0. For q < 0, the scale factor is 
always larger than in an EdS universe. 

For q < 0, the inflexion point indicating when the uni- 
verse starts decelerating is located at 



— = 2 1 / 3 |g| 1 / 3 . 
a 



(110) 



\/3-ln 



1 + V3 N 
V2 , 



1-9 



1/2 



(111) 



00, a c — > +00 and fV^V^oic — > 1.074. On the 



For q - r 

other hand for t —¥ —00, the scale factor behaves like 



'-ClnHnt 



"0 



(112) 




Fig. 10. Age of a BEC universe as a function of q > 0. The 
Einstcin-dc Sitter model corresponds to q = 0. 



Let us naively assume that all the content of the uni- 
verse is in the form of BEC dark matter and take fio = 1- 
The age of a BEC universe corresponds to the time to a t 
which a/ao = 1. For < q < 1, 



Hot = 



1 - 



<1 



l-q 



■ arctan 



1-9 



and for q < 0, 



H to = -U- 



1-9 



In 



(113) 



(114) 



-^0^0 — 77 • 



(115) 



As can be seen in Fig. [TQ1 the age of a BEC universe with 
q 7^ is smaller than the age of the EdS universe (q = 0). 

8.4. Contribution of radiation, baryons and dark energy 

Several observational results indicate that the universe is 
flat (k = 0) and that its present-day expansion is accelerat- 
ing (a > 0). Furthermore, the estimated age of the universe 
is to ~ 13.75 Gyrs which is about Hq 1 . If we assume that 
the universe is made only of a pressurclcss fluid and if we 
take A = 0, we are led to the EdS universe. However, this 
model leads to a decelerating expansion and predicts a too 
small age of the universe t = 2/3H ~ 9.3 Gyrs. If we as- 
sume that dark matter is a BEC, we again find that the 
present-day universe is decelerating. Furthermore, its pre- 
dicted age is even smaller than in the EdS universe (see Sec. 
18. 3[) . Therefore, we must invoke some form of dark energy 
to understand the present-day acceleration of the universe 
and its age. In the standard ACDM model, dark energy is 
due to the cosmological constant. On the other hand, to be 
complete, we must include the contribution of dark matter, 
baryonic matter and radiation. In that case, Eq. (|102j) is 
replaced by 



O 



B.D 



rad.O 



(a/a ) 3 (a/a ) 4 (a/a ) 3 - q 



(116) 

where £Ib,o, ^rad,o, ^dm,o and are the present-day val- 
ues of the density parameters of the baryonic matter, ra- 
diation, dark matter and dark energy. We have assumed 
that dark matter is in the form of BEC. The case of a 
pressureless dark matter is recovered for q = 0. Following 
Harko (2011), we have adopted the numerical values Hb.o = 
0.0456, n rad . Q = 8.24 10~ 5 , n DM ,o = 0.228 and Q. A = 0.726 
(Hinshaw et'al. 2009). 




Fig. 11. Scale factor a(t) as a function of time t for a uni- 
verse filled with baryonic matter, radiation, BEC dark mat- 
ter with k > 0, and dark energy. The different curves corre- 
spond to q = (pressureless dark matter, dashed), q = 0.01, 
q = 0.1, q = 0.303471 and q = 0.5. 
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2.5 r 




Fig. 12. Scale factor a(t) as a function of time t for a 
universe filled with baryonic matter, radiation, BEC dark 
matter with k < 0, and dark energy. The different curves 
correspond to q = (pressureless dark matter, dashed), 
q = -0.01, q = -0.1 and q = -0.5. 



The evolution of the scale factor a(t) is represented in 
Figs. QT] and [T^] for k > and k < respectively. In both 
cases, the universe is initially decelerating but finally en- 
ters in a phase of acceleration due to the effect of dark en- 
ergy (cosmological constant). In the case k > 0, previously 
treated by Harko (2011), the scale factor increases more 
rapidly than in the standard model (k = 0). Furthermore, 
the density of dark matter becomes infinite at a finite radius 
so that radiation does not dominate in the early universe. 
Accordingly, the universe starts from a(0) > at t = 0. In 
the case k < 0, the scale factor increases less rapidly than 
in the standard model (k = 0). Furthermore, the density of 
dark matter does not diverges when a — > so that radiation 
dominates in the early universe. Accordingly, the universe 
starts from a(0) = with an infinite (radiation) density. 

8.5. Theory of linear perturbations 

The theory of linear perturbations based on the Newtonian 
equations with pressure (|TT ]l -(|T9 |) has been performed by 
Reis (2003) in the general case, so that we shall directly 
apply his results to the present situation. Note that Harko 
(2011) considers other equations and obtains different re- 
sults. Nevertheless, our main conclusions will be the same. 
For an arbitrary equation of state, the evolution of the den- 
sity contrast is given by (Reis, 2003): 



5 - H [3(2i 



-eff 



2] 8 



-3H 2 



H C eff 



3w 2 
2 



Aw - - 
2 



3c 



2 2 

differ 



eff —A5 = 0, 



(117) 



where H = a/a is the Hubble constant, w = pb/pbC 2 , 
(Sp/Sp)/c 2 and c 2 = p' b {pb)/c 2 . These last quantities 



represent the velocity of sound normalized by the velocity 
of light. For a barotropic fluid, c 2 ^ = c 2 = p'(pt)/c 2 . On 
the other hand, from the Friedmann equations (|50")) - (|52"j) . it 
is easy to establish the relation 



1 + w 



= -3H(c 2 s - w). 



(118) 



We shall now focus on the equation of state ([55} . For this 



equation of state, we obtain w = kp b and c 
In that case, Eq. (|117[) reduces to 



.2 

eff 



6+2H(3w+l)5-'3H 2 



9w 2 



6w 



2wc 2 



2w. 



AS = 0,(119) 



where we have used w/H = — 3w(l + w) according to Eq. 
(|118p . Measuring the evolution in terms of a rather than 
in terms of t, and using the Friedmann equations (|8ip and 
(|82|) . we can rewrite Eq. (|119p in the form 



2 d 2 S 3, , d8 3, 2 



-12w+l)5- 



2wc 2 
H 2 a< 



Finally, according to Eq. (J5HJ), we have 
±1 



IV 



AS = 0.(120) 



(121) 



where the upper sign corresponds to k > and the lower 
sign corresponds to k < 0. If, for simplicity, we ignore the 
Laplacian term in Eq. (|120[) , we obtain 



d 2 S 
dR 2 



3 (1 + 3 ^JL_3 (9w 2 + 12w . 



1)6 = 0. 



For R- 
d 2 6 



R- 



dR? 



+oo, w — > and Eq. (|122[) reduces to 
2 dR 2 



(122) 



(123) 



This returns the usual equation (|55p of a pressureless fluid 
in the Einstein-de Sitter universe. Its solutions are given by 
Eq. (JSSD 

Let us now consider the case of small R. For k > 0, we 
have seen that the universe exists only for R > 1. Now, for 
R 1, w —> +oo so that the differential equation (|122|) 
is ill-defined. This means that we must start the perturba- 
tion analysis at t = U > 0. Harko (2011) remarks that the 
requirement c e // = c s < 1 (meaning that the velocity of 
sound must be smaller than the velocity of light) implies 
R > 3 1 / 3 and w < 1/2. Therefore, we shall assume that the 
universe starts at t* such that i?(i*) = 3 1 / 3 . According to 
Eq. dUl), we have 



3 
2 

For t 



KU = \/2 - arctan s/2 ~ 0.4589. 

t t , Eq. (|122p can be approximated by 

0. 



„d 2 5 I5„d5 111 
R 2 — - H R S 

4 dR 8 



dR 2 

and its solutions are 



SfxR- 



6otR- 



(124) 



(125) 



(126) 



The growth of the perturbations is faster than in a cold 
EdS universe (a 2 6 instead of a). On the other hand, for 
k < 0, the universe starts at R = for t — > — oo implying 
w = —1. In that case, Eq. (|122p can be approximated by 



R 



d 2 S 



dS 
~d~R 



dR 2 

and its solutions are 
5 cx i? 3 , 5 cx R. 



+ 3<5 = 0, 



(127) 



(128) 
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Contrary to a classical cold universe, the two solutions of 
the differential equation (|127[) are growing. Furthermore, 
the growth of the perturbations is faster than in an EdS 
universe (a 3 instead of a). Combining Eq. (I128[) with Eq. 
(JS5J), we get 



S(t) oc e 



3Kt 



6(t) oc e 



Kt 



(129) 



for t — > — oo. 

In conclusion, these simple estimates indicate that the 
growth of the perturbations is faster in a BEC universe than 
in a prcssurclcss universe. This conclusion was previously 
reached by Harko (2011) for k > based on different equa- 
tions. Our approach confirms this conclusion and extends 
it to k < (It may be recalled that the Laplacian term in 
Eq. (|119p has been neglected in our simple analysis. Now, 
we have seen in Sec. 16.21 that it is precisely this term that 
leads to an increase of the maximum growth rate in the 
(static) Jeans problem when k < 0. Therefore, the inclu- 
sion of this term in our analysis should enhance the growth 
of perturbations in the case of attractive self-interaction) . 

9. Dark fluid with a generalized equation of state 

9.1. Linear equation of state 

One possibility to understand the present-day acceleration 
of the expansion of the universe is to invoke a form of dark 
energy arising from a non-zero value of the cosmological 
constant A (see Sec. 18. 4[) . However, the physical meaning 
of A is not clearly understood. A major problem is that 
most quantum field theories predict a huge cosmological 
constant from the energy of the quantum vacuum, more 
than 100 order of magnitude too large. Therefore, other 
approaches have been developed to understand the phase 
of acceleration without invoking the cosmological constant. 
One possibility is to consider a "dark fluid" with a negative 
pressure. In this context, many workers (see the review of 
Peebles & Ratra 2003) have considered a linear equation of 
state 



p = ape , 



(130) 



with — 1 < a < 1. For this equation of state, the Friedmann 
equations ([80]). (j8l"j) and j82| with k = A = reduce to 



d P b i q^M , ^ 
— + 3-(l + a)p b 

at a 



0. 



a 
a 



4ttG 



(1 + 3a) p b , 



IttG 



-Pb- 



(131) 



(132) 



Equation (|131|) leads to the relation p b a 3 ^ 1+a ^ ~ 1. We note 
that the effect of a cosmological constant A is equivalent 
to a fluid with an equation of state p = —pc 2 (a = —1) 
since, in that case, p b is constant and can be written p b = 
Pa = A/8ttG. This equation of state leads to an exponential 
growth of the scale factor 



a e^ 



(133) 



On the other hand, for -1 < a < 1, Eqs. (|T3"Tj) and (fi"32l 
generate a model of the form 



ooct a/[3(i+o)] j 



a 3(1 + a)t ' 



(134) 



Pb 



6ttG(1 + a)H 2 ' 



(135) 



According to Eq. (|132I a). this universe is accelerating for 
— 1 < a < Q c = —1/3 and decelerating for a > a c . 
For a = a c , the scale factor increases linearly with time 
(a oc t). The EdS universe corresponds to a = 0. For 
the general model ()134[) - (|135p . the age of the universe is 
t Q = 2/[3(l + a)H ] where H„ = 2.273 10~ 18 s" 1 is the 
present-day value of the Hubble constant. We note that 
for the critical value a c = —1/3, the age of the universe 



is t 



13.95 Gyrs which is close to the value 



13.75 Gyrs predicted by the ACDM model. However, the 
deceleration parameter Q = —da/a 2 = (1 + 3a)/2 vanishes 
for a = —1/3 while its present-day value is ~ —0.5. We 
also note that the linear equation of state (|130p does not 
allow for a transition between a phase of deceleration and 
a phase of acceleration, while such a transition is likely in 
our universe. It may be therefore interesting to generalize 
this model. 



9.2. Generalized equation of state 

We consider a generalized equation of state of the form 

p= {ap + kp 2 )c 2 , (136) 

with — 1 < a < 1 and k positive or negative (the case a = 
— 1 is specifically treated in Appendix [U]) . In this approach, 
a single "dark fluid" combines the properties of a BEC dark 
matter described by the equation of state ([53")) and of a dark 
energy described by the equation of state (|130[) . We note 
that the equation of state of the BEC dark matter (oc p 2 ) 
dominates in the early universe where the density is high 
while the dark energy (oc p) dominates in the present-day 
universe where the density is low. The equation of state of 
the BEC also dominates at the scale of dark matter halos. 
This makes the study of this equation of state interesting. 
Another nice feature of this equation of state is that it 
admits fully analytical solution^. 

For the equation of state (|136[) . the Friedmann equation 
(|5D"|) becomes 



^+3%(l + a + fc P6 ) = 0. 
at a 

This equation can be integrated into 



Pb 



A(l + a) 



kA' 



(137) 



(138) 



where A > is a constant. For a — s- +oo, pi, ~ A(l + 
a)/a 3 ( 1+a \ When k > 0, the density exists only for a > 
a* = (kA) 1 /^ 1 ^. For a -> a», p b -> +oo. When k < 0, 
the density is denned for all a and has a finite value p b = 
(1 + a)/\k\ when a -> 0. 

Combining Eqs. ([5Tj) and (|136p . and taking A = 0, we 
obtain 



a 
a 



AnGp b 



(1 + 3a + 3kp b ). 



(139) 



12 Eq. ()136p can be written as pb = w(t)pbC 2 where w(t) is a 
function of time. Many authors have considered an equation of 
state of that form with some prescribed function w(t). In our 
model, w(t) is not an ad hoc function but it is explicitly given 
by wit) = a + kp b (t). 
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Let us define a critical density and a critical scale factor 
1 + 3a 



-2kA \ ^t+^T 
l + 3oJ 



(140) 



corresponding to a possible inflexion point (a = 0) in the 
curve a(t). When k > and a > —1/3, the universe is 
always decelerating (a < 0). When k > and a < —1/3, 
the universe is decelerating for pb > p c (i.e. a < a c ) and 
accelerating for pb < p c (i.e. a > a c ). When k < and 
a < —1/3, the universe is always accelerating (a > 0). 
When k < and a > —1/3, the universe is accelerating 
for pi, > p c (i.e. a < a c ) and decelerating for pb < p c (i.e. 
a > a c ). 

In order to determine the temporal evolution of a(t), we 
shall assume that the universe is flat (n = 0). Combining 
Eqs. (JHU) and (fHgj) . we get 



8nGA 



(1 



1/2 



Va 3 ( 1 +«) - kA 



(141) 



For a — » +oo, we recover the solution (|134l) - (|135l) . It is 
convenient to define a* = (|fc|^4) 1 /[ 3 ' 1+Q ^ and introduce 
R, = a/a*. In that case, the density is given by 



1 



1 



Pb 



\k\ i?3(l+") T l' 

and Eq. (|14ip can be rewritten 



R = 



KR 



Vi? 3 (!+«) T l' 



(142) 



(143) 



where the upper sign — corresponds to k > and the lower 
sign + corresponds to k < 0. On the other hand, we have 
defined the constant 



K = 



8ttG, 



3|fc 



■(1 



1/2 



For k > 0, the solution of Eq. (|143[) is 



3(1 



(144) 



-2ft, (145) 



^3(1+0) _ ! _ arctan y/R3(l+ a ) _ ! 



where the constant of integration has been set equal to zero. 
In this model, the universe starts at a finite time t = 
(see Sec. 19.41 for a revision of this statement) with a finite 
radius R(0) = 1 and an infinite density pb{0) = oo. For 
t -> 0, i? ~ 1 + [3/(4(1 + a))] 1 /3(ift)2/3. For i +oo, 
it asymptotically approaches the solution (|134l) - (ll35p i.e. 
i? - (3(1 + a)A"t/2) 2 /t 3 ( 1+a )l. when a > -1/3, the uni- 
verse is always decelerating (R < 0). When a < —1/3, 
the universe is decelerating for pb > p c (i.e. R < i?+ = 
(-2/(1 + 3a)) 1/[3(1+Q)1 ) and accelerating for p b < p c (i.e. 
R > R£). The time at which the universe starts accelerat- 
ing is 



3(1 + a) 



Kt K = 



-3(1 + a) 
1 + 3a 



— arctan 



-3(1 + a) 
1 + 3a 



(146) 



Some possible evolutions of R(t) corresponding to a < 
— 1/3, a = 1/3 and a > 1/3 are represented in Fig. [iUl 
For k < 0, the solution of Eq. (|143l) is 



ii 



k>0 


1 


/ 1 


- 


OK-1/3 y 








a = -1/3 






a>-l/3 

i 



Kt 

Fig. 13. Evolution of the scale factor in the case k > for 
different values of a (specifically a = —2/3, a = —1/3 and 
a = 1/2). 
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II 
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k<0 




/ i / 




,:< 1// / 








/ a = -1/3 






^^a>-l/3 


— 


i,i 





Kt 

Fig. 14. Evolution of the scale factor in the case k < for 
different values of a (specifically a = —2/3, a = —1/3 and 
a = 1/2). 



where the constant of integration has been set equal to 
zero. In this model, the universe starts from t — > —oo with a 
vanishing radius R(— oo) = and a finite density p&(— oo) = 
(1 + a)/|/c|. For t —> +oo, it asymptotically approaches 
the solution (flB^ -ipS ]) i.e. R - (3(1 + a)A"i/2) 2/[3(1+a)1 . 
When a < —1/3, the universe is always accelerating (a > 
0). When a > —1/3, the universe is accelerating for pb > p c 
(i.e. R< R~ = (2/(1 + Sa)) 1 /^ 1 ^) and decelerating for 
Pb < Pc (i-e. a > a c ). The time at which the universe starts 
decelerating is 



3(1 + a) 



Kt r 



3(1 + a) ^ I VT+3a+^+^ \ 



l + 3a 



V2 



For t — > — oo, R — >■ and the asymptotic expansion of Eq. 
(fT47)) yields 



i?(t) oc e 



Ki 



(149) 



^3(1+0) + i_i n 



1 + VW 11 



1 



i?3(l+a)/2 



3(1 + a ) (147) Some possible evolutions of R(t) corresponding to a < 
2 ' —1/3, a = 1/3 and a > 1/3 arc represented in Fig. ITU 
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9.3. Optimal parameters 

We can easily extend the analysis of Sec. l8.3l to the equation 
of state (|136[) . The density can be written 



Pb = Po 



where 



( JL)3(l+a) 

\ an ' 



kpo 



1 



kp 



implying kp 



g(l + «) 
l-<? 



(150) 



(151) 



We note that q is positive for k > and negative for k < 0. 
Furthermore, q = for k = 0, q — > 1 for k — > +oo and 
q —> — oo for k —> — (1 + a)/p . Using these results, the 
differential equation (|52"j) giving the evolution of the scale 
factor with the time can be rewritten 



a = HoV^o(l - q) 



vl/2. 



(±) 3{1+a} -q 



(152) 



We can make the connection with the previous notations 
by setting 



ao\q 



1/13(1+0)] 



K = H { 




(153) 



In that case, Eq. (|152p takes the form of Eq. (|143|) and its 
solutions are given by Eqs. (| 145|> and (|147p . Returning to 
the notations of this section, we find that the evolution of 
a{t) is given by Eqs. (|105[) - (|107[) where all the 3 are replaced 
by 3(1 + a). The other equations can be easily generalized. 

Let us assume that all the content of the universe is in 
the dark fluid with the generalized equation of state (|136[) 
so that f2rj = 1- Let us call to the time at which a /do = 1. 
For < q < 1, 



3(1 



for q < 0, 



Hntn — 



3(1 



1 



1 - 



■ arctan 



1 



(154) 



and for q = 0, 



3(1 



(155) 



(156) 



These expressions basically give the age of the universe to 
as a function of q and a. 

The deceleration parameter is defined by 



da 1 + 3w(t) 



(157) 



where w(t) = Pb/pbC 2 . For the equation of state (|130D . 
w{t) = a + kpi,. Using Eq. (|151|) . we find that the present- 
day value of the deceleration parameter can be expressed 
in our model as 



to = 



3a + 2q+l 
2(1^9) ' 



(158) 



The observed value of Qo is close to —0.5. If we take Qo = 
-1/2, Eq. ([155)1 yields a = -(g + 2)/3. On the other hand, 



in the standard model, the age of the universe is about 
13.75 Gyrs. If we take to = H^ 1 (which is close to this 
value) and use the previous relation between a and q, Eqs. 
(|154j) and (|155[) determine the optimal value of q, then a. 
For < q < 1, the optimal value of q is given by 



l + q 



* .arctan. Z 1 "* 



1-9 



0. 



(159) 



We find q = 0.303471 and a = -0.76782367. For q < 0, the 
optimal value of q is given by 



1 




q + V^q) = o. 



(160) 



We find q = -0.383589 and a = -0.53880367. The evo- 
lution of the scale factor a(t) in a universe filled with a 
dark fluid with equation of state (|136[) is represented in Fig. 
1151 for the optimal values of (q, a) obtained previously It 
gives a relatively good agreement with the standard ACDM 
model based on a non-vanishing value of the cosmological 
constant. 




Fig. 15. The full line represents the scale factor a(t) as 
a function of time t for a universe filled with a dark 
fluid with equation of state (|136|) with optimal parame- 
ters (q = 0.303471, a = -0.76782367) and (q = -0.383589, 
a = —0.53880367). The long-dashed line corresponds to 
the standard model (see Figs . [TT1 and [T2"|) . By construction, 
these three curves give the same values of the age of the 
universe to ~ ^o" 1 an( l °f the present-day deceleration pa- 
rameter Qo ~ —0.5. The model with k > is decelerating 
then accelerating. Using Eqs. (|146[) and (|153[) . we find that 
the acceleration starts at Hot c ~ 0.2 which is earlier than 
the value H n t c ~ 0.5 predicted by the standard model. 
The model with k < is always accelerating. We have 
also indicated the model corresponding to a linear equa- 
tion of state (<7 = 0) with a = —1/3 yielding a/ao = Hot. 
It gives the correct age of the universe but has a vanish- 
ing deceleration parameter. Finally, we have indicated the 
(pressurelcss) EdS model corresponding to a = q = yield- 
ing a/a = (3H t/2) 2 / 3 . It does not give the correct age 
of the universe and is decelerating instead of accelerating. 
The other models with a = and g / do not do bet- 
ter. This shows that we must combine the two equations 
of state ap and kp 2 to obtain a good agreement with the 
standard model, i.e. we need an equation of state with two 
parameters. 
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9.4. Theory of linear perturbations 

In this section, we extend to the case of the generalized 
equation of state (|136[) the theory of linear perturbations 
developed in Sec. 18. 51 For this equation of state, w = a+kpi, 
and c^jj: = c 2 = 2w — a. According to Eq. (|142[) . we have 



w = a ± 



1 + a 



R 3(l+ a ) zp 1 ' 



(161) 



where the upper sign corresponds to k > and the lower 
sign corresponds to k < 0. The generalization of Eq. (|122l) 
is 



if 



d 2 S 
dR 2 



dS 
dR 



-(9w 2 + 12w + 1 - 8a - 6a 2 )S = 0. 



(162) 



For R —t +oo, w — >• a and Eq. (|162[) reduces to 
d 2 S 3, s „dS 3 



R 1 



dR 2 



(l-a)R—--(a + l)(3a + l)S = 0. 



(163) 



This is the equation for the density contrast corresponding 
to the linear equation of state (|13Q[) . Its solutions are 



S ocR 



l+3a 



5<xR-^ 1+a K 



(164) 



Let us now consider the case of small R. For k > 0, 
the requirement c e // = c s < 1 implies R > [(3 + a)/(l — 
Q! )] 1 /[3(i+«)] anc j w < (a + l)/2. Therefore, we shall assume 
that the universe starts at t* such that i?(i») = [(3+a)/(l— 
a )]i/[3(i+a)]_ According to Eq. ([IS5]) . we have 



3(1 + a) -KU 



2(1 + a) 2(1 + a) , , 

' arctanj \ ' . (165) 



2 ~"~ V I- a "V 1-a 

For t — > i„, Eq. (|162[) can be approximated by 

and its solutions are 



(166) 



-ll + 15n±\/-135« 2 -90a + 1009 

5± oc s 



(167) 



For k < 0, the universe starts at i? = for t — oo imply- 
ing w = — 1. In that case, Eq. (|162p can be approximated 

by 

R2 4l^ - 3 (! + 2a ) R T7i + 3 (! + + 3c 0<5 = °, ( 168 ) 
dit^ ait 

and its solutions are 

SocR 3{1+a \ 5<xR 1+3a . (169) 
10. Conclusion 

Following the proposal of Bohmer & Harko (2007) and oth- 
ers (see a short historic in Chavanis 2011), we have assumed 
that the dark matter in the universe is a self-gravitating 
BEC with short-range interactions, and we have theoreti- 
cally explored the consequences of this hypothesis. For the 



sake of generality, we have considered the case of positive 
and negative scattering lengths. 

At the level of dark matter halos, a positive scattering 
length, equivalent to a repulsive self-interaction generat- 
ing a positive pressure, is able to stabilize the halos with 
respect to gravitational collapse. This leads to dark mat- 
ter halos without density cusps, equivalent to polytropes 
of index n = 1 (more generally, the barotropic equation 
of state is fixed by the form of the self-interaction so that 
other configurations are possible) . Alternatively, if the scat- 
tering length is negative, equivalent to an attractive self- 
interaction generating a negative pressure, the dark matter 
is very unstable and collapses above a very small critical 
mass M max = 1.0127i/ y/\a s \Gm (Chavanis 2011). When 
these ideas are applied to an infinite homogeneous cosmic 
fluid (Jeans problem), it is found that a negative scatter- 
ing length can increase the maximum growth rate of the 
instability and accelerate the formation of structures. The 
virtues of these results could be combined by assuming that 
the scattering length changes sign in the course of the evo- 
lution. It could be initially negative to help with the for- 
mation of structures and become positive (due to a change 
of density, magnetic fields, radiation,...) to prevent com- 
plete gravitational collapse. The mechanism of this change 
of sign is, however, unknown so that this idea remains 
highly speculative. However, some terrestrial experiments 
have demonstrated that certain atoms can have negative 
scattering lengths, that their scattering length can depend 
on the magnetic field, and that it is possible in principle 
to manipulate the value and the sign of a s (Fedichev et al. 
1996). Therefore, a BEC is a serious candidate for which 
the pressure can be positive and/or negative (see Appendix 
iBl for further remarks about the values of the BEC param- 
eters). 

At the cosmological level, we have constructed mod- 
els of universe composed of BEC dark matter with at- 
tractive or repulsive self-interaction. We have first studied 
the academic situation where the universe is made only of 
BEC dark matter. A BEC universe with positive scatter- 
ing length, having a positive pressure, is not qualitatively 
very different from a classical Einstein-de Sitter universe. 
It also emerges at a primordial time t = from a big- 
bang singularity where the density is infinite, and under- 
goes a decelerating expansion asymptotically equivalent to 
the EdS universe. A difference, however, is that the ini- 
tial scale factor a(0) is finite. On the other hand, a BEC 
universe with negative scattering length, having a negative 
pressure, markedly differs from previous models. It starts 
from t — > — oo with a vanishing radius and a finite density, 
has an initial accelerating expansion then decelerates and 
asymptotically behaves like the EdS universe. This model 
universe exists for any time in the past and there is no 
big-bang singularity. When we add the effect of radiation, 
baryonic matter and dark energy (via the cosmological con- 
stant), the picture is different. In that BEC universe 
with attractive or repulsive self-interaction starts from a 
singularity at t = where the density is infinite. It first 
experiences a phase of decelerating expansion followed by 
a phase of accelerating expansion. For k — > we recover 
the standard ACDM model but for k ^ 0, the evolution 
of the scale factor in a BEC universe can be substantially 
different. The model with k > expands more rapidly than 
the standard model (fc = 0). The initial scale factor is fi- 
nite (a(0) > 0) and the radiation never dominates. The 
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model with k < expands less rapidly than the standard 
model. The initial scale factor vanishes (a(0) = 0) and the 
radiation dominates leading to a decelerating expansion. 
In both models, the dark energy dominates at large times 
leading to an accelerating expansion. Finally, we have con- 
sidered a "dark fluid" with generalized equation of state 
p = (ap + kp 2 )c 2 having a component p = kp 2 c 2 similar 
to a BEC dark matter and a component p = ape 2 mimick- 
ing the effect of the cosmological constant (dark energy). 
We have found optimal parameters (a, k) that give a good 
agreement with the standard model. We have studied the 
growth of perturbations in these different models and con- 
firmed the previous observation of Harko (2011) that the 
density contrast increases more rapidly in a BEC universe 
than in the standard model. 

In conclusion, the idea that dark matter could be a BEC 
is fascinating and probably deserves further research. 

Appendix A: Solution of the differential equation 

tfSHD 

In this Appendix, we provide some analytical solutions of 
the differential equation (|58p . If we define 



3fc 4 
2fi£ 



A = 



3fc^ 
2k 2 



this equation can be written 



d?l_2_df_ 
da 2 2a da 



1 ffi 

,2 



7^ + ^7 l- + \a 2a 



4-3 7 



1 / = 0. 



With the change of function f(a) = a 3 ^ 4 g(a), we get 



d 2 g 1 dg 1 



da 2 



da 



a 16 , 



Assuming a^O and defining x = a a , we obtain 



>d 2 g , dg 1 



dx 2 



dx 



r X/a 



Xx 2 



25 
16 



9 = 0. 



(A.3) 



(A.4) 



We have not been able to find the general solution of this 
differential equation in terms of simple functions. We there- 
fore consider particular cases. We first consider the case 
p = 0. Then 



dx 2 



X ^- + A- (xx 2 - — 
dx a 2 \ 16 



0. 



(A.5) 



For A > 0, defining y = \[Xxja andp = 5 /4a, the foregoing 
equation becomes 



2 d 2 g , _dg 



+ yJ + (y 2 -P 2 )g = o. 



dy 2 dy 



(A.6) 



This is a Bessel equation whose regular solutions are g 
J± p (y). Coming back to the original variables, we get 



f(a) cx a 3/4 J ± ^ 



(A.7) 



For A < 0, defining y 
becomes 



i*9 



2 u_y + dg_ 
dy 2 dy 



f^Xx/a and p = 5/4a, Eq. (|A.5j) 
{y 2 +P 2 )g = 0- (A.8) 



This is a Bessel equation whose regular solutions are g 
I± p {y). Coming back to the original variables, we get 



f{a) « a 3 / 4 I ± 5 



-Aa" 



a 



(A.9) 



For a = and fx = 0, Eq. (|A.2[) becomes 

d 2 f 1 df 1 „ 
dai-2-a£ + ^^-^ f = - 

Looking for solutions of the form f(a) cx a n we find that 



(A.10) 



/(n)otuW, (A.ll) 
Finally, for a = and p ^ 0, Eq. (|A.2[) can be rewritten 



fl 1 df , 1 (l t 
da 2 



-T + 4(^ + a-i)/ = o. 

a da a z \a J 



(A.12) 



Proceeding as before, we find that the solution of this equa- 
tion is 



f(a) cx a 3/4 J ±i 



±i\/25-16A 



(A.13) 



Appendix B: Estimate of the BEC parameters 

In this Appendix, we derive some constraints on the BEC 
parameters. The parameter k appearing in the equation 
of state of a BEC is defined by k = 2ira s h 2 /m 3 c 2 



(A.l) 



(A 2) an d the present-day critical density is (p c )o = 3Hq/8ttG. 
Introducing the dimensionless parameter A = a s mc/h 
(Chavanis 2011), we obtain the general relations 



HPc)o = \{H t p ) 2 (^) 



eV/c 2 



2.45710 



-10 



and 
A = 5.09110" 



fm eV/c 2 



(B.l) 



(B.2) 



On the other hand, the radius of a BEC dark matter halo 
is given by (Arbey et al. 2003, Bohmer & Harko 2007, 
Chavanis 2011): 



R = 7T 

so that 

( — 

Vkpc 



a s % 
Gm 3 



2 \ 1/2 



5.951 10 4 



A^\ 1/2 : 
Gc 



cV/c 2 



(B.3) 



(B.4) 



Therefore, the typical size of the dark matter halos deter- 
mines the ratio A/to 4 . Taking R = lOkpc, we obtain 
(to/(cV/c 2 ))/A 1 / 4 = 4.94. Now, we remark that the pa- 
rameter k(p c )o, given by Eq. (|B.1|) . only depends on A/to 4 . 
Using the previous estimate, we get k(p c ) = 4.126 10~ 13 . 
We conclude that the dimensionless parameter q, defined 
by Eq. is very small. 

This estimate, which does not rely on any free param- 
eter, seems to indicate that the pressure of the BEC dark 
matter is totally negligible at the scale of the cosmic fluid 
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(since q ~ like in a pressureless universe) while it is impor- 
tant at the scale of galactic halos (since R 2 ~ A/m 4 ). This 
seems to be bad news for the BEC cosmology. In fact, this 
result implies either that (i) the BEC dark matter can be 
treated as a pressureless fluid, like the ordinary dark mat- 
ter, at the cosmological scale (i.e. in the Friedmann equa- 
tions) or that (ii) the BEC parameters (in particular the 
scattering length a s ) are different in the homogeneous cos- 
mic fluid (before the Jeans instability) and in the dark mat- 
ter halos (after the Jeans instability) . This observation may 
be a hint that the value of the scattering length changes in 
the course of time and that the scattering length of the 
bosons in the dark halos is not the same as in the cosmic 
fluid (because the density is different). This corroborates 
our remark (Chavanis 2011) that the sign of the scattering 
length may change too. Clearly, the determination of the 
BEC parameters, or the constraints that they must satisfy, 
is certainly a very important step for the validation of the 
BEC dark matter hypothesis. 



Appendix C: Equation of state p = (—p + kp 2 )c 2 



References 



In this Appendix, we consider the equation of state (| 136[) 
with a = — 1 and k ^ 0. In that case, Eq. (|137[) can be 
integrated into 



Pb 



3k ln(a/a„) 



(C.l) 



where a* is a constant. Physical solutions require that k > 
and a > a*. Setting R = a/a*, the Friedmann equation 
with k = A = can be written 



R= 2 -K* 
3 V^R 



(C.2) 



where K = (2nG/k)^ 2 . The solution of Eq. ((021) is 

R(t) = e^ V3 . (C.3) 



For t -> 0, R ~ 1 + (Kt) 2 ^. The curve R(t) presents an 
inflexion point at R c = y/e, (pb) c = 2/ (3k) and Kt c = 
(1/2) 3 / 2 . The expansion is decelerating for t < t c and ac- 
celerating for t > t c . 

If we introduce the present-day values po and ao of the 
density and scale factor, the density can be rewritten 



Pb = 



Po 



3gln(a/arj) 



1 



(C.4) 



where q = kpo. We have a*/ao = exp(— l/3g) and K = 
(3fl /iq) 1/2 H . Using Eq. (fU3j) the evolution of the scale 
factor is given by 



a(t) 
ao 



1 , 3S"i ^1/3, it ,\2/3 



(C.5) 



The time at which the universe starts accelerating is Hot c = 
(q/GQo) 1 / 2 ■ The age of the universe, corresponding to the 
time t = to at which a/ao = 1 is Hgto = 2/(9q) (we have 
taken fio = 1). Finally, the present-day value of the decel- 
eration parameter (|157[) is Qq = (3q — 2)/2. 
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